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ABSTRACT 

We describe a new series of band-limited image masks for coronagraphy that 
are insensitive to pointing errors and other low-spatial-frequency optical aberra- 
tions. For a modest cost in throughput, these "eighth-order" masks would allow 
the Terrestrial Planet Finder Coronagraph (TPF-C) to operate with a pointing 
accuracy no better than that of the Hubble Space Telescope. We also provide 
eighth-order notch filter masks that offer the same robustness to pointing er- 
rors combined with more manageable construction tolerances: binary masks and 
sampled graded masks with moderate optical density requirements. 

Subject headings: astrobiology — circumstellar matter — instrumentation: adap- 
tive optics — planetary systems 
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1. INTRODUCTION 

Coronagraphy holds great promise for imaging extrasolar planetary systems, even ex- 
trasolar terrestrial planets only ~ 10~ 10 times as bright as their host stars (e.g. Kuchner & 
Spergel 2003a). However, finding extrasolar terrestrial planets at contrast levels of ~ 10~ 10 
using any of the present image mask designs requires either pointing accuracies at the level 
of a fraction of a milliarcsecond (Kuchner & Traub 2002; Kuchner & Spergel 2003b) or 
apodization in the pupil plane (e.g. Aime et al. 2002; Kasdin et al. 2003), which generally 
carries a high penalty in throughput and inner working angle (but see also Guyon 2003; 
Traub & Vanderbei 2003). We offer a new series of band-limited image masks that can pro- 
vide high contrast levels without pupil apodization because they are intrinsically insensitive 
to pointing errors and other low-order aberrations. We also provide notch filter versions of 
these masks that may be easier to build to the necessary tolerances. 



2. EIGHTH-ORDER MASKS 

2.1. Band-limited Masks and Notch Filter Masks 

Here we summarize the basic definitions of band-limited masks and notch filter masks 
stated by Kuchner & Traub (2002); Kuchner & Spergel (2003b). We will focus on linear 
masks, described by functions of a single variable, x. One-dimensional band-limited and 
notch filter masks can be combined to create a wide variety of two-dimensional masks. 

An ideal linear image mask can be described by a function, M(x), called the mask 
function. In our simple model of the interaction between masks and light, the mask function, 
also called the mask's amplitude transmissivity, multiplies the electric field phasor of the 
incoming beam. The intensity transmissivity of a mask, |M(a;)| 2 , multiplies the intensity of 
the beam. We will also refer often to the mask function, the mask intensity transmissivity, 
and to the Fourier transform of the mask function, 

M{u) = J M(x)e~ 27riux dx. (1) 

Kuchner & Traub (2002) showed that if M(x) is a notch filter function, i.e., 

M(u) = for e/2 < \u\ < 1 - e/2, (2) 
where e sets the undersizing of the Lyot stop, and if 

/ M(u) du = 0, (3) 

J-e/2 
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then the mask denned by M(x) will completely remove all on-axis light in an ideal corona- 
graph with a uniform entrance pupil. Kuchner (2005) showed that notch filter masks are the 
only trivially achromatic masks that completely remove on-axis light in a one-dimensional 
or separable two-dimensional coronagraph. Masks we can construct without amplifying the 
beam or manipulating its phase are necessarily limited to < M(x) < 1. A band-limited 
mask is a notch filter mask with M(u) = for \u\ > e/2. 

We aim to find notch filter mask functions, M(x), that provide deep suppression of 
light near the optical axis, not just at the optical axis. We will first derive new band- 
limited masks and then follow the recipes in Kuchner & Spergel (2003b) to generate useful 
notch-filter masks based on them. 



2.2. Blocking Slightly-Off-Axis Light 

Understanding the off-axis behavior of an ideal coronagraph with a band-limited mask 
is easy. A coronagraph with a band-limited image mask attenuates the intensity of an image 
of a point source located at an angle i by a factor of |M(a;)| 2 compared to the image the 
source would have if the image mask were removed while the Lyot stop remained in place 
(see the Appendix). In an ideal coronagraph with a band-limited mask, the point spread 
function (PSF) is independent of the position of the source with respect to the optical axis; 
only the attenuation varies with x. Hence, we can describe the way a band-limited mask 
attenuates sources near the optical axis, including the target star, simply by expanding M(x) 
about x — 0. 

If the first important term in this expansion is quadratic in x, the intensity attenua- 
tion will vary as x 4 . Borrowing the language of interferometry, we might say such a mask 
produces a fourth-order null. For a demonstration of why this interferometric terminol- 
ogy is appropriate, consider the nulling coronagraph described by Levine et al. (2003), which 
monochromatically synthesizes a particular band-limited mask with a fourth-order null using 
beam combiners. 

All of the band-limited mask designs and notch filter mask designs illustrated in Kuchner 
& Traub (2002), Kuchner & Spergel (2003b), and Kuchner (2005) have fourth-order nulls. 
For example, all of the popular 1 — sine™ family of masks are fourth order; 1 — sinc n kix/n 
(l/6n)(A;ia;) 2 . But we can design band-limited masks and notch filter masks with nulls of 
any order, (3, by the methods described below, if (5 is a multiple of 4. 

The order of the null dictates the sensitivity of the mask to optical aberrations that 
effectively spread the light from a target source around some region of the sky near the 



-4- 



optical axis. Tip-tilt error (caused by pointing error, for example) is the simplest low-order 
aberration for us to model and a term that can easily dominate a coronagraph design's error 
budget. A pointing error of Ad will cause an intensity leak proportional to (A6) 13 . A mask 
that is insensitive to pointing error will also defeat other low-order aberrations like defocus, 
coma and astigmatism to some degree, though some low-order Zernike terms contain mid- 
spatial-frequency tails that may leak through. Mid-spatial-frequency errors are problematic 
for any coronagraph design because by definition they coincide with the search area; no mask 
or stop can block them without also blocking light from the planet. Shaklan & Green (2005) 
discuss the effects of low-order aberrations in a coronagraph with an eighth-order mask in 
detail. 

The fractional leakage through a mis-pointed coronagraph with a band-limited mask is 
simply 

J JI(x + A6,y)\M(x,y)\ 2 dxdy 

ffl(x,y)dxdy ' 1 } 

where I(x, y) is the source intensity, i.e., the stellar disk, and AO is the instantaneous pointing 
error. For a fourth-order linear mask, the instantaneous fractional intensity leakage is 

T _ fl 4 + 48^ 2 (Ag) 2 + 128(Afl) 4 

where is the angular diameter of the star and 9iw is the inner working angle of the 
mask, defined by \M(6 IW )\ 2 = 1/2. To derive this expression, we made the approximation 
that M(x) = x 4 ; we have corrected a numerical error in Equation 17 of Kuchner & Spergel 
(2003b). If we assume A6 is distributed in a Gaussian with standard deviation UAe, and 
o"A6» >> 9*, then we find that the mean leakage is 

(L) = 1.5(a A0 /e IW ) 4 . (6) 

So if we assume that we can tolerate a leakage of (L) < 3 x 10~ 8 , and that AO is much larger 
than the angular radius of the star, we find that we must center the star on the mask to an 
accuracy of 

a Ae < 0.012 6 IW . (7) 

Though it is easiest to interpret in terms of pointing error, this Gaussian blurring can also 
serve as a crude model of the effects of other low-order aberrations. 

For an eighth-order mask approximated as M(x) = x 8 , the instantaneous fractional 
intensity leakage is 

_ 7 6 8 + 1120 el(A6) 2 + 17920 fl 4 (Afl) 4 + 57344 ^(A^) 6 + 32768(Afl) 8 

65536 6\ /2 ' ( } 
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the corresponding mean fractional leakage is 

(L) = 52.5(a Ae /6 1/2 ) s , (9) 

and the pointing requirement for leakage (L) < 3 x 1CT 8 is 

a Ae < 0.070 IW , (10) 

a factor of ~ 6 improvement over the a A e tolerance for fourth-order masks. A coronagraph 
designed to find extrasolar terrestrial planets like the Terrestrial Planet Finder Coronagraph 
(TPF-C) might need 9 IW = 60 milliarcseconds (mas). This requirement implies a pointing 
tolerance of a A o < 0.72 mas using a fourth-order mask or a A Q < 4.2 mas using an eighth- 
order mask. For comparison, the Hubble Space Telescope points to cr A g ~ 3 mas (Burrows 
et al. 1991). 

Eighth-order masks can also provide high-contrast images of extended sources, though 
relaxing the pointing tolerance depletes some of this power. For a fourth-order mask, Equa- 
tion 5 shows that the extent of a central source begins to matter when > (8/3) (A^) and 
the cross term begins to dominate. For an eighth-order mask, Equation 8 shows that the 
extent of a central source begins to be important when 6* > (4/7) (A^). In the TPF-C exam- 
ple above, these limits correspond to = 1.9 mas for a fourth-order mask and 6* = 2.4 mas 
for an eighth-order mask; a solar-type star at 10 pc is about 1 mas in diameter. In other 
words, a TPF-C design with an eighth-order mask may be slightly better suited for the 
closest target stars than one using with a fourth-order mask even with its relaxed pointing 
tolerance, depending on the wings of the actual distribution of pointing errors. 



3. CONSTRUCTING THE MASKS 

To design an eighth-order band-limited mask, we can create a linear combination of two 
fourth-order band-limited masks weighted so that the term responsible for the quadratic leak 
cancels; i.e., 

d 2 



dx A 



0. (11) 



x=0 



For example, we can add a term of the form C(l — cos(/c2^)), otherwise known as a sin 2 mask, 
to any 1 — sine™ mask to create a new mask with d 2 M(x)/dx 2 \ x= o, while still satisfying 
Equation 3. If we start with a mask of the form 1 — sinc n kix/n « (kix) 2 /(6n) and add 
C(l — cos(/c 2 a;)) ~ (C/2)(A; 2 a;) 2 , we find that to produce an eighth order mask, we require 
that C = -(1/3^(^/^2)2. 
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However, we do not want to add a sin 2 mask of just any random spatial frequency. We 
would prefer a frequency within the bandwidth of the original mask so that we don't suffer 
an undue throughput penalty; i.e., hi needs to be < h\. In order to minimize |C|, we should 
pick a frequency at exactly the edge of the band; i.e., k 2 — k±. With this constraint, we find 
C= -l/(3n). 

Of course, adding mask functions can violate the requirement that M(x) < 1. To ensure 
M(x) < 1, we can renormahze the mask by multiplying M(x) by a constant, N, equal to 
the inverse of its maximum value. 

Putting everything together and using physical units yields a series of eighth-order 
band-limited masks, 



Mbl{x) = N 



3n — 1 . nxe ( 1 \ nxt 
— sine — + — cos ■ 



(12) 



nXmaxf \3n J ^rnaxf _ 

where / is the focal ratio at the mask and X max is the longest wavelength at which the mask 
is to operate. Figure 1 shows M(x) for the first few linear masks in the series. The n = 3 
design offers a good compromise between the large sidelobes of the n — 1 mask and the 
higher inner working angle-bandwidth product of the n — 5. 

The ringing in these image masks reduces their effective throughputs. The amplitude 
of the additional ringing introduced by the cosine term in Equation 12 falls off slowly with 
n, so simply increasing n does not help much. 

Fortunately, we can create another series of eighth-order masks with less ringing by 
combining two 1 — sine™ masks instead of a 1 — sine™ mask and sin 2 mask using the same 
procedure we used to construct Equation 12: 



M BL (x) = N 



I — m . , irxe m . ttxs 
— sine — + — smc 



(13) 



I lX max f I 7Tl\ ma xf _ 

This series of masks has less ringing than the series described by Equation 12. It is 
parametrized by two integer exponents, / and m; we assume I > m. Figure 2 shows M{x) 
for m — 1 and I =2-5. The m = 1 and / =2-3 masks have throughput similar to the n = 3 
cosine mask. Using large values of m and / reduces the ringing further, but it also reduces 
the Lyot stop throughput. 

Figure 3 compares the intensity transmissivity, |M(x)| 2 , for the 1 — sine 2 fourth-order 
mask and the m — 1, I — 3 eighth-order mask. While the 1 — sine 2 mask has an inner 
working angle of 9iw = (1.448/e)(A/D), the m — 1, I — 3 eighth-order mask has an inner 
working angle of 9iw = (1.788/e)(A/D). The m — 1, I — 3 mask offers a good compromise 
between ringing and throughput, and also reaches 100% transmissivity at its first maximum, 
a critical region for planet searching; we recommend this mask for TPF-C. 
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distance from center of mask (X I e D) 



Fig. 1. — Eighth-order band-limited mask functions described by Equation 12 for n = 1 — 5. 



- 8- 




distance from center of mask (UeD) 



Fig. 2. — Eighth-order band-limited mask functions described by Equation 13 for m = 1, 
/ = 2-5. 
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distance from optical axis (k I D) 




distance from center of mask {X I e D) 



Fig. 3. — Intensity transmissivities for the 1 — sine 2 fourth-order mask, the n = 3 eighth-order 
mask, and the m — 1, I — 3 eighth-order mask. Coronagraph throughput and distance from 
optical axis were calculated with e = 0.6. The m — 1, I — 3 eighth-order mask, recommended 
for TPF-C, has 100% transmissivity at its first maximum. 
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Consider a TPF-C design with 9iw = SX/D using a linear fourth-order mask. This 
coronagraph has a bandwidth of e = 0.4 and a nominal Lyot stop throughput of 1 — e = 
0.6 (see Kuchner & Spergel 2003b). This fourth-order design probably requires some mild 
apodization of the Lyot stop to ameliorate leakage due to low-order optical aberrations, 
reducing the throughput to 0.35. Keeping 6iw = 3X/D but switching to a linear m — 1, 
/ = 3 eighth-order mask would mean working at a bandwidth of e = 0.596 and a Lyot stop 
throughput of 1 — e = 0.404. Coronagraphs with eighth-order masks should not require any 
Lyot stop apodization. 

In other words, our analysis suggests that eighth-order masks combined with un-apodized 
Lyot stops perform about as well as fourth-order masks combined with apodized Lyot stops 
in terms of throughput and robustness to pointing errors. An alternative way to provide 
robustness to pointing errors is to use a shaped-pupil coronagraph (Kasdin et al. 2003; Van- 
derbei et al. 2003a,b). But the throughput offered by an eighth-order linear mask is still 
better than the typical throughput of a shaped pupil coronagraph at a given inner working 
angle, particularly when a shaped pupil coronagraph is used with a hard-edged image mask, 
which increases its effective inner working angle (Kuchner 2005). 



4. EIGHTH-ORDER NOTCH FILTER MASKS 

The functions described by Equations 12 and 13 can be used in a variety of ways, e.g., 
to make linear masks (M(x,y) = Mbl(x)), radial masks (M(r) = Mbl(j')), or separable 
masks (M(x,y) = Mbl{x)Mbl(v))- However, all band-limited masks are necessarily smooth 
graded masks. Notch-filter masks offer even more design freedom and need not necessarily be 
smooth, making them potentially easier to manufacture than band-limited masks (Kuchner 
& Spergel 2003b). 

Notch filter masks affect starlight and planet-light the same way as band-limited masks; 
only their low-spatial frequency parts contribute to starlight suppression. Consequently, in 
an eighth-order notch filter mask, only the low-frequency part needs to satisfy Equation 11. 
In other words, 



d 2 

—— (the low-frequency part of Mix) 
dx l 



6/2 

{-27riu) 2 M{u) du = 0. (14) 

x=0 J-e/2 



Equivalently, we can say that an eighth-order notch filter mask satisfies Equations 2, 3, and 
also 



re/2 

/ u 2 M{u)du = 0. (15) 

J-e/2 
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To find masks that meet these criteria, we can use a technique similar to the one 
employed in §3. Since any linear combination of fourth-order notch filter mask functions will 
automatically satisfy Equations 2 and 3, we will start by writing 

(16) 

where M no tch 4A i x ) and M no tch 4B (%) represent different fourth-order notch filter mask functions 
and N ensures that M n otch a ( x ) < 1- To construct a notch filter mask that exhibits eighth- 
order behavior, we need to weight the linear combination so that the new notch filter function 
also satisfies Equation 15. In other words we will find the constant C by substituting 
Equation 16 into 15. 

_ P-e/1 u2M notch 4A (U) du 

— • (17) 

J_ e/2 u 2 M notchiB (u) du 
This constant should be negative; it should also satisfy \C\ < 1 

By combining fourth-order notch filter functions and using the solutions to Equation 17, 
we can construct a variety of eighth-order notch filter masks analogous to the variety of 
eighth-order band-limited masks. For example, we can make a family of eighth-order notch 
filter masks using the 1 — cos and 1 — sine™ fourth-order notch filter functions. We can also 
design low-ringing eighth-order notch filter masks using the 1 — sine" 1 and 1 — sine' notch 
filter functions. To be consistent with §3, we will refer to the various eighth-order notch filter 
masks by the exponents of their constituent functions (n, m, I, ... etc.). In the following, we 
provide example calculations for making eighth-order binary and graded notch filter masks 
using the m — 1, I — 3 design. 



M notchs (x) = N M notch4A (x) + CM notchiB (x) 



4.1. Eighth-order Binary Masks 

Notch filter masks can be designed to be binary: everywhere either completely opaque 
or completely transparent. A simple way to make such a binary mask is to assemble a mask 
from a collection of identical parallel stripes, where any arbitrary notch filter mask function 
provides the width of each stripe. In other words, each stripe is defined by 

jCf (rr ,a _ J 1 where V < M notch (x) \ min f 

M stripe (x,y) - ^ Q elgewhere _ (18) 
and the mask function is 

oo 

Mf)i nar y(^X ', ?/) ^ ^ M stripe 

(X, y-jX m inf), (19) 

j=—oo 
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where \ m i n is the shortest wavelength in the band of interest. 



If we like, we can use the band-limited mask functions described by Equations 12 or 13 
in place of M notch (x), resulting in a mask formed of continuous curves. However, sampled 
binary masks may prove to be easier to manufacture since their features are not as small 
near the optical axis. We will construct here an eighth-order sampled binary mask. Such 
a mask can be made entirely from rectangles of opaque material. Debes et al. (2004) have 
demonstrated the construction of sampled fourth-order masks using e-beam lithography. 

Fourth-order sampled masks are defined by the following prescription (Kuchner & 
Spergel 2003b): 

M notch4 (x) = M samp4 (x) - M , (20) 

where 



M, 



samp4 



M, 



sampi 



P(x) * (^Mbl 4 (x) Ax S(x - (k + QAx^j , 
P(u) ( M BL4 (u) * 6 ( u - k/Ax)e~ 2mu < Al 



(21a) 
(21b) 



and 



eD/(2A) 



/oo 
M BLi (u)P(u) du = M BLi (x) *P(x) . (22) 

Here, M B l a represents any fouth-order band-limited mask function, k ranges over all integers, 
and * indicates convolution. The sampling points are offset from the mask center by a fraction 
( of Ax. The kernel, P{x), can represent the "beam" of a nanofabrication tool. It should be 
normalized so that P(x) dx = 1, and P(x) must be everywhere < l/(Ax), so M samp4 (x) 
remains < 1. The constant M ensures that the mask satisfies Equation 3. Though the 
sampled mask is derived from Mbl 4 (x), the function being sampled is M B l 4 (x) — M . 



Combining Equation 16 and Equation 20, we have 

M notchs (x) = N M notch4A (x) + CM notch 

A' {M samp4A (x) -M 0a ) + C(M samp4B (x) - M 0b ) 



(23) 
(24) 



where M, 



sampan 



[x] are sampled versions of the fourth-order band-limited functions M B l 



4A,B 



x 



described by Equation 21b. The constants M 0ab and C ensure that M notc h 4A B (x) satisfy 
both Equation 3 and Equation 15, The constant C is given by 



C = 



J e J e 2 /2 u 2 P(u)M BL4A (u) du 
p_!* 2 v?P{u)M BL4B {u) du 



(25) 
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Figure 4 shows a plot of the function N[(M BUa (x) - M 0a ) + C(M BLiB (x) - M 0s )] to 
illustrate how ( may be chosen. This example uses the m — 1, I — 3 sampled eighth-order 
mask, meaning M B l 4A (x) = 1 — sinc(:r) and Mbl 4B (x) = 1 — sincere). To guarantee that 



will choose ( = Co, to make the central rectangles contiguous. 

The bandwidth of a mask should be chosen conservatively; e.g., A max should be some- 
what larger then the longest wavelength where the detector is sensitive so a filter with a 
finite slope can remove all the extraneous light. Band-limited masks and notch filter masks 
leak light at wavelengths longer than X max ] notch filter masks also leak light at wavelengths 
shorter than \ rn in- At a fixed inner working angle, increasing \ max necessitates increasing 
e, and thereby decreasing the throughput. Decreasing \ min means spacing the stripes and 
samples in a notch filter mask closer together. 

For the m — 1, I — 3 mask with Oiw = 3A ma;c /-D, spacing Ax = X m i n f, and bandpass 
0.5-0.8 jum, we find that M 0a = 0.00630889, M 0b = 0.01882618, C = -0.33935486, and 
Co = 0.25941279. Table 1 lists normalization constants and sampled mask parameters for 
eighth-order notch filter masks of various inner-working-angles using a top-hat kernel, P(x) = 
(D I \ min )U(xD I \ m in) i an d 0.5-0.8 /iin bandpass. 

If the resolution of our nanotool is ~ 20 nm, we require a telescope with an f/115 or 
slower beam (see Kuchner & Spergel 2003b). The physical size of an entire mask is generally 
a few hundred diffraction widths. A 1" x 1" mask would consist of > 440 vertically repeating 
segments, where each segment is < \ m inf = 57.5 /im wide. This coronagraph design would 
have a Lyot stop throughput of 40%. Figure 5 shows an example of what an m — 1, I — 3 
linear eighth-order binary mask would look like. 

A similar mask can also be made by replacing Equation 18 with 



as shown in Kuchner & Spergel (2003b). In this design, the notch-filter function is reflected 
vertically for each segment. Manufacturing a mask like the one shown in Figure 5 could 
substantially reduce the writing time for e-beam lithography and opportunities for write 
errors. 



M notc h(x) > 0, the parameter ( must be in the range \(\ < Co, where Co is defined by the 
condition M BL4A (( \ min f) + CM B L 4B (Co^minf) = M 0a + CM 0b . For our binary mask, we 





elsewhere. 



(26) 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 4. — An example of the function N[(M BLiA (x) - M 0a ) + C(M BL4B (x) - M 0b )} for an 
m — 1, I — 3 sampled mask. Choosing C = Co allows us to create a binary mask of contiguous 
stripes. Choosing ( = allows us to create a graded mask with the most favorable optical 
density requirement. 



Table I. SAMPLED EIGHTH-ORDER MASK PARAMETERS 



n 


N a 


&IW {^max/ D) 


e 




Mo B 


C 


Co 


D b 

. max 


1 


0.966115405054 


3 


0.453 


0.01092315 


0.01631996 


-0.67562610 


0.25945500 


8.004 




0.960497496515 


4 


0.340 


0.00616927 


0.00923329 


-0.67167174 


0.25954657 


9.012 




0.959860814806 


5 


0.272 


0.00395309 


0.00592117 


-0.66985758 


0.25958912 


9.791 


2 


0.999927046667 


3 


0.487 


0.00632078 


0.01883265 


-0.34113344 


0.25944228 


7.969 




0.999967637078 


4 


0.366 


0.00357716 


0.01068997 


-0.33769399 


0.25953913 


8.969 




0.999991487843 


5 


0.292 


0.00227903 


0.00682024 


-0.33609563 


0.25958482 


9.757 


3 


0.994355716928 


3 


0.533 


0.00505061 


0.02250740 


-0.22964671 


0.25941276 


7.860 




0.992967898001 


4 


0.400 


0.00284944 


0.01275232 


-0.22635050 


0.25952294 


8.868 




0.992249764357 


5 


0.320 


0.00182510 


0.00818418 


-0.22484879 


0.25957404 


9.649 


4 


0.999920502046 


3 


0.578 


0.00445606 


0.02640459 


-0.17384571 


0.25937869 


7.744 




0.999959029620 


4 


0.434 


0.00251629 


0.01499184 


-0.17065228 


0.25950368 


8.751 




1.000006649135 


5 


0.347 


0.00160975 


0.00961517 


-0.16919646 


0.25956186 


9.534 


/ i 


^for m — 1) 
















2 


1.865785172445 


3 


0.557 


0.00825681 


0.01646433 


-0.50505400 


0.25942680 


7.923 




1.862096989484 


4 


0.412 


0.00452972 


0.00904463 


-0.50274761 


0.25953456 


8.977 




1.856230853161 


5 


0.334 


0.00298028 


0.00595415 


-0.50180101 


0.25957920 


9.710 


3 


1.434216871605 


3 C 


0.596 


0.00630889 


0.01882618 


-0.33935486 


0.25941279 


7.882 




1.429552473250 


4 


0.447 


0.00355642 


0.01063737 


-0.33669458 


0.25952307 


8.890 




1.427349701514 


5 


0.357 


0.00227076 


0.00679929 


-0.33546973 


0.25957440 


9.674 


4 


1.312506672966 


3 


0.637 


0.00540801 


0.02147409 


-0.25623997 


0.25623997 


7.813 




1.308947497039 


4 


0.478 


0.00305104 


0.01215389 


-0.25348152 


0.25951079 


8.819 




1.306220598720 


5 


0.382 


0.00195033 


0.00778076 


-0.25221408 
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Fig. 5. — Simulated high magnification picture of an m — 1, 1 — 3 linear eighth-order sampled 
binary mask. Dark areas are completely opaque and white areas perfectly transmissive. See 
Debes et al. (2004) for optical microscope photographs of an actual linear fourth-order 1-sinc 2 
binary mask. 



a Normalization constant for C = Co an d f^min sampling. 
b For a graded mask with ( — 0. 
"Suggested for TPF-C. 
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4.2. Eighth-order Graded Masks 

Smooth graded band-limited image masks have produced suppression of on-axis monochro- 
matic light at the level of a few times 10 -9 in the laboratory (Trauger et al. 2004). However, 
construction errors probably still limit the broad-band performance of these masks. We 
suggest that sampled graded masks may be easier to construct than smooth graded masks. 
Graded masks produce large phase errors, but it may be possible to correct the phase of these 
sampled masks using transparent strips of varying thickness. Also, as Kuchner & Spergel 
(2003b) pointed out, sampled masks can be designed so that, unlike smooth masks, they do 
not require their darkest regions to be perfectly opaque. This flexibility limits the demands 
on the lithography tool used to make the masks. The 1 — sine 2 mask with 6jw = 2.9 \ max /D, 
e = 0.4, can be built with a maximum optical density of 4. The sin 2 mask with e = 0.4 can 
be built with a maximum optical density of 3. 

When we design eighth-order graded notch filter masks, we can reduce the required 
maximum optical density by beginning the sampling at ( = 0, so long as the spacing between 
the samples is large enough to straddle the valleys shown in Figure 4. Choosing Ax = \ m inf 
satisfies this condition for all of the masks listed in Table 1. Figure 6 shows a graded 
version of the m — 1, I — 3 eighth-order mask described in §4.1. The mask is defined by 
M(x, y) = M notch (x); its optical density is — log 10 \M notch (x)\ 2 . To make the darkest stripe of 
the mask as transparent as possible, we chose ( — 0. With this choice, the darkest stripe of 
the mask has optical density —2 log 10 \M notc h(0)\ ~ 7.882. Table 1 lists the maximum optical 
densities (O.D. max ) of sampled graded masks with £ = 0. 




Fig. 6. — Simulated low and high magnification pictures of an m — 1, 1 — 3 eighth-order sam- 
pled graded mask with ( = 0. The low magnification picture (left) shows ~ 400 diffraction 
widths; each stripe has uniform shading and is f\ m in wide. 
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5. SUMMARY 

We offered a series of eighth-order masks that are relatively insensitive to tip-tilt errors 
and other low-spatial-frequency aberrations; in a coronagraph using one of these masks, the 
r.m.s. pointing error only needs to be managed to a few milliarcseconds, no better than 
the pointing accuracy of the Hubble Space Telescope. Eighth-order notch filter masks retain 
most benefits of using fourth-order masks — broadband capabilities, reasonably high through- 
put, and small inner working angle — permitting extremely-high dynamic range coronagraphy 
suitable for terrestrial planet finding using a popular optical layout. 

In particular, we suggested a binary mask designed for TPF-C at 0.5-0.8 /im composed 
of opaque strips whose shapes are described by Equation 24 with m — 1, I — 3, e = 0.596, 
N = 1.434216871605, M 0a = 0.00630889, M 0s = 0.01882618, C = -0.33935486, and Co = 
0.25941279. This mask provides 40% Lyot stop throughput and requires an f/115 or slower 
beam, assuming the mask can be manufactured with an r.m.s. accuracy of 20 nm. The 
r.m.s. pointing accuracy required for achieving starlight suppression of 10 -10 with this mask 
in the search area is a^e ~ 4.2 milliarcseconds for stars of diameter up to ~ 2.4 mas. If the 
mask is used on a telescope with better pointing accuracy, it can achieve contrast levels of 
10~ 10 on targets with even larger diameters. 

We also provided a graded version of this design, whose amplitude transmissivity is 
described by Equation 24 using the above parameters but with ( = 0. This mask offers 
the same performance as the above binary version, but it allows easier e-beam fabrication 
because it only requires optical densities < 7.882. Other eighth-order masks can provide less 
ringing at the cost of inner working angle or Lyot stop throughput. 

We thank Stuart Shaklan and Joseph Green for helpful conversations and for delaying 
the publication of their paper on low-order aberrations in coronagraphs with eighth-order 
masks until this paper was ready. M. J.K. acknowledges the support of the Hubble Fellowship 
Program of the Space Telescope Science Institute. J.C. and J.G. acknowledge support by 
NASA with grants NAG5-12115, NAG5-11427, NSF with grants AST-0138235 and AST- 
0243090, the UCF-UF Space Research Initiative program, and the JPL TPF program. 

A. APPENDIX 

We will prove for a monochromatic coronagraph with a notch filter mask, a binary 
entrance aperture of finite size, and a Lyot stop that is perfectly opaque everywhere the 
entrance aperture is opaque, that 1) the PSF shape is the absolute square of the Fourier 
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transform of the Lyot stop amplitude transmissivity independent of the position of the source 
on the sky, and 2) the PSF is attenuated by the intensity transmissivity of the band-limited 
part of the mask evaluated at the source position. Kuchner & Traub (2002) demonstrated 
this principle for a sin 2 mask; this more general proof applies to any two-dimensional notch 
filter mask. 

As usual, we will examine a coronagraph comprising an entrance aperture, A, an image 
mask, M, and a Lyot stop, L, each of which is represented by a complex-valued function. 
We will use the notational conventions of Kuchner & Traub (2002) and Kuchner & Spergel 
(2003b): letters with hats represent image plane quantities. The image-plane coordinates 
are x = (x,y) and the pupil-plane coordinates are u = (u,v). 

Monochromatic light propagates through the coronagraph as follows. 

1) An incoming wave incident on the entrance aperture creates a field with amplitude E(u). 
When an incoming wave interacts with a stop or mask, the function representing the mask 
multiplies the wave's complex amplitude. So after the wave interacts with the entrance 
aperture, the amplitude becomes A(u) ■ E(u). 

2) After the entrance aperture, the beam propagates to an image plane, where the new field 
amplitude is the Fourier transform of the pupil plane field amplitude, A(x) *i?(x); * denotes 
convolution. In this plane, the beam interacts with the image mask, and the field amplitude 
becomes M(x) • (i(x) * E(x)). 

3) Next, the beam propagates to a second pupil plane, where the field amplitude is M(u) * 
(A(u) ■ E(u)). In this second pupil plane, the wave interacts with a Lyot stop, changing the 
field amplitude to F(u) = L(u) • [M(u) * (A(u) ■ E(u))}. 

4) At last, the beam propagates to the final image plane, where the final image field is 
F(x), the Fourier transform of F(u). For a point source, the intensity of the final image is 
proportional to the absolute value of this quantity squared. 

The final image field, -F(x), and its Fourier transform are linear functions of A(u), L(u), 
and also M(u). This last property allows us to study masks by decomposing them into 
Fourier components, computing F(u) or F(x) for each one, and then summing the final field 
amplitudes back together. 

Consider a point source providing a field -E'(x) = 5(x — xi) in the plane of the sky 
and a harmonic mask function M(u) = S(u — Ui). The field after the entrance pupil 
is A(u) exp(— 2ni u ■ Xn), and the field in the first image plane is A(x — x x ). The field 
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after the image mask is exp(27riux • x)A(x — xi). The field in the second pupil plane is 
A(u — u x ) exp(— 27ri(u — u x ) • x x ). The field after the Lyot stop is 

F(u) = L(u)A(u - u 1 )e- 27Ti(u - Ul) - Xl for a harmonic mask. (Al) 

Let A be binary (everywhere equal to 1 or 0) and let A represent the support of A 
and £ represent the support of L. If C G A, then there is some set V G 1Z 2 for which 
L(u)A(u — ux) = L(u) for u x G V. If A is finite in extent, then there is also some set 
QeTZ 2 for which L(u)A(u - u x ) = for u x G Q. 

Under these circumstances, there are three kinds of harmonic image masks: 

Ux G V: For these harmonic masks, the field after the Lyot stop is uniform in amplitude 
with a phase gradient x x . 

ux G Q: For these masks, the field inside the Lyot stop is zero. 

Ux G" (PUQ): For these masks, the field after the Lyot stop does not have uniform amplitude. 

These three kinds of harmonic masks correspond to the three kinds of virtual pupils illus- 
trated in Figure 6 of Kuchner & Traub (2002). 

A band-limited mask is defined to be a continuous sum of harmonic masks of the first 
variety; 

M(x)=/ M(ux)e 27riuix dui. (A2) 
JuieP 

A notch filter mask is defined to be a continuous sum of harmonic masks of the first and 
second varieties; 

M(x)= / M(ui)e 2mui - X dui. (A3) 

iuie(PuQ) 

Combining this expansion and Equation Al using the linear property of F(u) described 
above, we find that in a coronagraph with a notch filter mask, the field amplitude after the 
Lyot stop is 

F(u) = / M(ui)L(u)A(u - ui)e- 2 " (u - Ul) x Mui (A4) 

iuiG(PUQ) 

= I M(ui)L(u)e- 2m(u - Ul) ' Xl dm (A5) 

= L(u)e- 2 ™ X1 I M(ux)e 27riui xi dm. (A6) 
iuieP 

(A7) 
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To interpret this equation, let us define the band-limited part of M(x) as 

M BL (x)=/ M( Ul )e 27riuiX d Ul . (A8) 

Now we can write 

F(u) = M BL (x 1 )L(u)e- 2 ™' Xl for a notch filter mask. (A9) 

The final image field is the Fourier transform of this quantity, F(x) = M^£,(x 1 )I/(x — xi), 
and the final image intensity is the absolute square of the Fourier transform of this quantity, 

|F(x)| 2 = |M BL (x 1 )| 2 |L(x - xi)| 2 for a notch filter mask. (A10) 

In other words, for a notch filter mask, the PSF shape is |L(x)| 2 , independent of Xi, the 
position of the source on the sky. The PSF is attenuated by a factor |Mbl(xi)| 2 , the ampli- 
tude transmissivity of the band-limited part of the mask evaluated at the source position. 
The band-limited part of a notch filter mask can generally be found by applying a low-pass 
filter to the mask function. 
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